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We derive new positivity bounds on spin-dependent parton 
distributions in multicolored QCD. They are stronger than 
Soffer inequality. We check that the new inequalities are sta- 
ble under one-loop DGLAP evolution to higher normalization 
points. 

1. In spite of the vast phenomenological data pro- 
viding information about the quark and gluon distribu- 
tions in nucleons we are still far from the complete knowl- 
edge of all twist-two distributions, especially of the spin- 
dependent ones. Therefore various positivity bounds on 
parton distributions are very helpful. In particular, an 
important role is played by Soffer inequality Q that con- 
straints the transversity quark distribution which is not 
measured yet. At first sight the set of already known 
positivity bounds for twist-two distributions is complete 
and no enhancement can be derived on general grounds 
without solving the dynamics of QCD. However, if one 
takes QCD in the limit of large number of colors N c then 
it is possible to obtain new stronger inequalities. 

It is well known that in the large N c limit g] the 
baryons are described by mean field (Hartree) equations 
H , which can be written in terms of bilocal meson fields 
with all possible spin and isospin quantum numbers. Al- 
though the corresponding effective bilocal meson action 
is not known, a number of results have been obtained for 
the large N c baryons without specifying the dynamics. 
Usually these results use the spin-flavor symmetry 
of the solution of the Hartree equation describing the 
nucleon and A resonance: this solution is invariant un- 
der simultaneous space, spin and isospin rotations (but 
not invariant under a pure flavor rotation). Owing to 
this spin-flavor symmetry the baryon states (described 
by this mean field solution) appear as "rotational exci- 
tations" with spin S and isospin T connected by con- 
straint S = T In particular, this approach allows 
to describe nucleon (S = T = 1/2) and A resonance 
(S = T = 3/2). 

This picture of the large N c baryons is extensively 
used in the Skyrme model and various chiral quark- 
soliton models |g] . But some general model- independent 
results have been extracted for the large N c baryons us- 
ing only the spin-flavor symmetry ||. This includes, for 
example, the identities for the mass splitting of baryons 
and the large N c classification of spin and isospin struc- 
tures of various form factors. In this paper we study 
quark distribution functions relying only on the large N c 
limit and on the spin-flavor symmetry of the large N c 
mean field solution. 



We first present our main result — the new positivity 
bounds which can be considered as a serious enhance- 
ment of the Soffer inequality and then give the derivation 
and check the stability of the bounds under one-loop the 
DGLAP evolution. 

2. Let us first formulate the new positivity bound. 
At finite number of QCD colors the unpolarized (<?/), 
longitudinal polarized (A^qf) and transverse polarized 
(At5/) quark (antiquark) distributions are known to be 
constrained by the following set of inequalities 



\ a lQ/\ < 9/ , 
I A T ?/ 1 < 9/ , 
q f + A L q f > 2|A T <7/| 



(finite N c ) 



(1) 



The last inequality is known as Soffer inequality [EJ. 
These inequalities hold for each quark flavor /. In the 
leading order of the large N c limit we have for the u and 
d distribution functions of the proton Jl(]|ll|] : 

q u = q d , A L q u = -A L q d , A T g„ = -A T q d . (2) 

Therefore the set of usual inequalities ([[]) takes the fol- 
lowing form in the large N c limit 



\A L q u \ < q u , 

\A T q u \ < qu , 
q u - \A L q u \ > 2\A T q u 



(large A^ c old bounds) (3) 



Note that in the last inequality (||) we have the absolute 
value of A^q u in contrast to Soffer inequality at finite 
N c (|l|). Actually we combined two Soffer inequalities for 
f = u,d (Q) and we made use of the large A^ c relation 
A L q u = -A L q d (see eq. (§)). 

In this paper we shall derive the following enhancement 
of inequalities @ 

\q u + A L q u > 2\A T q u \ , 



(large iV c new bounds) (4) 



3<7„ - A L q u > . 
Let us introduce compact notations 
A L q u A L q d 



L = 



T = 



A T q u 



Qd 
A T q d 



O 
O 



qu qd 
Then our new result ([|) becomes 
h + L>2\T\, 



(new) 



(5) 



(6) 



L > 0, 



1 



whereas the old inequalities (g) in the limit N c — > oo 
result in 



m<i, 

ILI + 2ITI < 1 



(old) 



(7) 



Graphically the old and new regions of the allowed values 
in the L, T plane is shown in Fig. [I]. 




FIG. 1. Softer inequality at large N c (the rhombus hatched 
with vertical lines) versus new large N c inequalities (the tri- 
angle hatched horizontally) in terms of variables L and T (see 
eq. (§). 

3. Now we turn to the derivation of our main re- 
sult (H|). In order to consider the parton distributions 
in the large N c limit we start from the quark correlator 
(B2\ip(x2)'ip{xi)\Bi) . The Hartree picture of the large N c 
nucleon described in || leads to the following form of this 
correlation function 

(Bi , Pi | ^+ fl (x? , xi )V> S2 f 2 (z° ) x 2 ) I B 2 , p 2 } = 

N c J d 3 Xe l( - P2 ~ Pl) * J dRD* Bi (R)D B2 (R) (8) 



x i 



x®, xi 



X,x 2 -X ){R- V 



Here F S2 f 2tSl f 1 (x\ — x^^x,^) is the quark correlation 
function (ip + (xi)ip(x2)) in the background field of the 
solution to the Hartree equation, Si are Dirac spinor in- 
dices and / are SU (2) isospin indices. The x° — y° time 
dependence expresses the fact that we deal with a static 
solution. Solution F S2 f 2 _ Sl f 1 (x\ — x%, xi, x 2 ) also violates 
the SU (2) flavor and the 3-space translation invariance. 
Therefore in eq. (||) we used plane waves e lp X and the 
rotator wave functions Db(R) depending on SU(2) ma- 
trix R in order to construct baryon states with given spin, 
isospin and momentum. 

Since the exact form of the effective Hartree hamilto- 
nian corresponding to the large N c limit is not known we 
do not know the function F S2 f 2tSl f 1 {x\ — :Kj,Xi ) X2) ei- 
ther. However, following || we assume that this solution 
has the spin- flavor symmetry: 



SU(2) S 

pill 



C SU(2) &a 



SO(3, 1)lo* 



(9) 



The transformations of ( 2 ) spin— flavor perform a simul- 
taneous rotation both in the isospin and in the usual 3D 
space. The invariance of the mean field solution under 
SU (2) spin _ flavor means that 



S S2S ' 2 {R)Rf 2 f. 2 {R 1 )f[hS s ' lSl { R X ) 
x F 4fiAfl (fi ~ x °2, 0(R)X!, 0(i?)x 2 
= -^2/2,81/1(^1 -x§,xi,x 2 ). 



(10) 



Here we use notation S S2S i 2 (i?) for spin rotations and 
0{R) for space rotations associated with isospin matrix 
R. 

The quark distribution functions can be expressed 
through the correlation function (ph in the standard way. 
Using the spin-flavor symmetry C P-Oj ) and performing the 
integration over R in (^) we easily arrives at the large 
N c relations (||). Moreover, using expression (||) one can 
show that any quark parton distribution function fi(x) 
can be represented in the large A^ c limit as follows 



fi(x) = Sp[Tip(x)] . 



(11) 



Here I\ and p{x) are matrices with SU{2) flavor indices 
and Dirac spin indices so that we deal with the trace of 
8x8 matrices (4 sp i„ x 2j SOS pj n = 8). Matrix p(x) can be 
expressed through function F S2 f 2tSl f 1 (n) as follows: 



p{x) 



Nc 
2tt 



dx° / d 3 xi / d 3 x 2 



^ 3 



(27T) 



P e ip(xi 



0„,3 



x e 



- i A.M w -/ ) l + 7V jr(x ;Xi[X2) l + f7 



-x 2 ) 



(12) 



Matrix p(x) is determined by the dynamics of the large 
N c QCD and is not known. But p(x) is universal for all 
twist- two parton distributions. 

Matrices Ti in eq. ( pd[ ) are determined by the specihe 
type of the distribution function fi(x). Due to the large 
iV c constraints (0) we can restrict to three independent 
distribution functions fi (x) listed in the table below with 
the corresponding 8x8 matrices Ti 



i 


fi 







Qu + Qd 


1 


L 


A L q u - A L q d 




T 


A T q u - A T q d 





Although matrix p(x) can be found only if one solves 
the large N c QCD, one can establish its general proper- 
ties: 
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1) p is a hermitean matrix 

P + = P; 

positive in the matrix sense 



(13) 



p>0, (14) 

2) p lives in the subspace of the projector (1 +7°7 3 )/2 

P7°7 3 =7°7 3 P = P, (15) 

3) p commutes with i^j 2 + r 3 due to the spin-flavor 
symmetry (|To|) 



[p, ( l7 V + r 3 )] = 0. 



(16) 



One could have an impression that our statements vi- 
olate the Lorentz and isotopic invariance. But actually 
everything is correct: the mean field solution for the large 
N c nucleon has the mixed spin-flavor invariance (^J) but 
due to the choice of the third axis for the boost of the 
infinite momentum frame we are left only with the axial 
spin-isospin rotations. As a result our eqs. (fi"5|), ( |l6| ) are 
invariant only with respect to the U(l) axial spin-isospin 
rotations. However, we stress that the final results for 
the parton distributions are certainly both Lorentz and 
isospin invariant. 

Using the properties (PI), ( |l6"l) of matrix p(x) and dis- 
crete P, C, T symmetries one can obtain the following 
representation for it 



1 + 7°7 3 



[ Cl l + c 27 5 r 3 + c 37 5 (7 1 T 1 + 1 2 t 2 )} (17) 



7 5 t 3 , 



7 5 7 1 t\ 



with some coefficients Cj. Matrices 7 °7 3 
7 5 7 2 t 2 appearing in this equation commute with each 
other and have eigenvalues ±1. Therefore they can 
be diagonalized simultaneously. Restricting the consid- 
eration to the 4-dimcnsional subspace of the projector 
(1 + 7°7 3 )/2 (see (|l5|)) we can diagonalize p as follows 

p = diag(ci - c 2 + 2c 3 , c\ + c 2 , C\ + c 2 , c% - c 2 - 2c 3 ) . 



The requirement of positivity (J14 
constraints on coefficients Cj 

ci - c 2 > 2|c 3 | 



ci + c 2 > . 



(18) 

leads to the following 



(19) 



(20) 



Inserting the explicit representation for p ( |l7j) into the 
general expression for parton distributions (|11[) we im- 
mediately express the parton distributions fi(x) through 
the coefficients Cj, and the constraints (19), (EOT) on these 
coefficients together with (||) lead to bounds (Q) on par- 
ton distributions. 



4. Now we would like to check that our new positivity 
bounds are preserved by the one-loop evolution to higher 
normalization points p with DGLAP kernels taken in the 
leading order of the large iV c limit. The derivation follows 
the ideas of papers |l^,|l3| . 

In the leading order of the large N c limit the one-loop 
evolution of the quark distributions is not affected by the 
gluon distribution and the evolution equations for quark 
distributions listed in the above table take the form 



P-j^fo = P\\ ® fo , 

M^/i = P|| ® /i , 

P-irfT = Pr <S> h , 
dp 



(21) 



where the singlet unpolarized quark distribution fo and 
the longitudinally polarized isovector distribution fi, 
have the same one-loop DGLAP kernel Pi | . 

We want to show that if for some numbers a, (3 at a 
given normalization point the following inequalities are 
satisfied 



«/l < fo , af L + (3f T < f c 



(22) 



then the evolution to higher normalization points pre- 
serves these inequalities. Indeed, introducing functions 

fa = fo ~ afh , fa,/3 = fo - otf L ~ (3f T , (23) 
we can rewrite the evolution equations as follows 



P 



dp \ fa 



U,p \ _( Pt (P\\-Pt) 







Pi 



fa,p 



Here all kernels are positive |14| : 

P T > 0, Pi I -P T > 0, Pi > 



(24) 



(25) 



so that the positivity of functions P a ,/3 and F a is pre- 
served during the evolution to larger p (note that the 
subtraction terms in the evolution equations do not vi- 
olate this positivity This completes the proof of 
the evolution stability of inequalities (^2|) for arbitrary 
coefficients a, (3. In particular, we can choose for a, (3 
the coefficients appearing in our bounds (||). Thus the 
evolution stability of inequalities (P is established. 

5. Comparing our new bounds with the phenomeno- 
logical data one should keep in mind that the allowed 
regions shown in Fig. [I] correspond to the leading order 
of the large N c expansion and that the 1/N C corrections 
can cause the phenomenological distribution functions to 
exceed the "triangle boundary" of Fig. 0. Indeed, if we 
take the GRV [|l5| parametrization for q u + qa and GRSV 
|16| for Al<7 u — A^qd then we shall see that our inequality 



3 



\A L q u - A L q d \ 1 

: < r (large N c ) 

q u + q<i 3 

is essentially modified by the 1/N C corrections: 



\A L q u - A L q d \ 
q u + qd 



0.6. (GRSV) 



(26) 



(27) 



Actually this deviation from our large N c bounds agrees 
with the fact that in chiral quark soliton models the 1/N C 
corrections to gA are known to be large |l7| ]. 

The large N c inequalities derived in this paper also 
hold for antiquark distributions. Since our knowledge 
about A^qu — A^qd and Axq u — Axqd is still scare the 
new bounds can be helpful in pinning down the polar- 
ized antiquark distributions (note that model estimates 
indicate that 1/N C corrections to polarized antiquark dis- 
tributions are rather small Jll| ) . The knowledge of polar- 
ized antiquark distributions is important for the analysis 
of data on semi-inclusive and Drell-Yan processes with 
polarized nucleons. 

Apart from pure phenomenological applications we 
think that the new large N c bounds can be used as a 
consistency test for various models of polarized parton 
distributions. 
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